function [ S ] = OutputProc(M,P,L)
%OUTPUTPROC Return to spatial domain
S = ifft2(ifftshift(M).*exp(1i*P),'symmetric');
S= uint8(S);
%image adjustment to values 0 - 255
% S = imadjust(S,[single(min(min(S)))/255 ; single(max(max(S)))/255] , [single(min(min(L)))/255 ; single(max(max(L)))/255] );